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Abstract 

00 . 

The vacuum energy density of electromagnetic field inside a perfectly con- 
ducting wedge is calculated by making use of the local zeta function technique. 



(N 

CN , This regularization completely eliminates divergent expressions in the course 

of calculations and gives rise to a finite expression for the energy density 
in question without any subtractions. Employment of the Hertz potentials 
for constructing the general solution to the Maxwell equations results in a 
considerable simplification of the calculations. Transition to the global zeta 
function is carried out by introducing a cutoff nearby the cusp at the origin. 
Proceeding from this the heat kernel coefficients are calculated and the high 
^ ■ temperature asymptotics of the Helmholtz free energy and of the torque of 

the Casimir forces are found. The wedge singularity gives rise to a specific 
high temperature behaviour ~ T 2 of the quantities under consideration. The 
obtained results are directly applicable to the free energy of a scalar massless 
field and electromagnetic field on the background of a cosmic string. 
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I. INTRODUCTION 



The Casimir calculations are usually accomplished by two methods: a global approach, 
when from the very beginning the total vacuum energy of quantum fields are calculated 
for given boundary conditions, and a local consideration, when first the space density of 
the vacuum energy is found. The both methods are in agreement if the integration of 
the vacuum energy density gives the total Casimir energy. As known, the first example 
where this condition is not satisfied was a perfectly conducting wedge, i.e., two perfectly 
conducting plates that are crossed at given angle. The calculation of the Casimir energy for 
these boundary conditions has the following peculiarities. The spectrum of electromagnetic 
oscillations inside the wedge is continues and it does not depend explicitly on the opening 
angle of the wedge. The boundaries are flat therefore there is no problem with additional 
divergences due to the boundary curvature. However, on the cross line of the plates the 
smoothness of the boundary is obviously violated due to the cusp. Till now a correct 
definition of the total vacuum energy of the fields inside the wedge is lacking. It has been 
calculated only the spatial density of vacuum energy inside the wedge. It turns out that this 
density possesses a non-integrable singularity ~ r -4 in cylindrical coordinates r, 8, z with 
the z axes along the cross line of the plates. For the first time these calculations have been 
done probably in Ref. by making use of a simple method of images and later in Ref. |2| 
by applying the Green's function method. In recent publications 0,(H in the framework of 
the Schwinger source theory it was shown how to avoid some mathematical inconsistency 
taking place in previous calculations @J5| . 

The geometry of space analogous to that between the crossing plates also encounters 
in other physical problems, in particular, in the cosmic string theory (the space outside 
the string) and under consideration of quantum fields at finite temperature on the Rindler 
space-time. 

The geometry of space in the background of a cosmic string is determined by the solution 
of the Einstein equations with the energy-momentum tensor of matter defined by the string. 
When the string is located along the z axes and possesses the linear mass density /z, then the 
space outside the string is isomorphic to the manifold C^-^ x R 1 , where M 1 is an infinite line 
along the z axes and instead of a plane perpendicular to this axes one has two dimensional 
cone C2n-(p with the angle deficiency <fi = 8TcGfi, G being the gravitational constant. The 
cone of an angle f3 becomes a plane when (3 = 2ir. Certainly the boundary conditions for 
fields inside the wedge W a x R 1 and on the space Cp x IR 1 are in the general case different. 
For the cone Cp the periodic boundary conditions are natural. For fields inside the wedge 
the boundary conditions are determined by the relevant physical origin of the fields and by 
the physical characteristics of the boundary. 

The manifold Cp x IR 2 arises under consideration of fields on the Rindler space-time 
at finite temperature j3~ x . In order to calculate the Helmholtz free energy and entropy of 
fields on this space the local zeta function technique was proposed ||. A nontrivial problem 
encountered here is the construction of required analytic continuation of this zeta function. 
Such a continuation has been accomplished separately for lower and higher eigenvalues of 
the Laplace operator in the problem at hand. 

As far as we know the Casimir energy for a perfectly conducting wedge at finite tem- 
perature was considered only in Ref. [[jj by the multiple scattering method. In that paper 
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unusual temperature dependence of the Casimir forces in this problem was noted. These 
forces, tending to reduce the angle a between the plates, diminish when taking into account 
the first temperature correction. 

The present article seeks to show that the most simple and consistent (from the math- 
ematical point of view) approach to the calculation of the Casimir effect for a perfectly 
conducting wedge is to use the local zeta function technique. In order to construct such 
spectral function one should find not only the spectrum of the regarding operator but also 
its eigenfunctions. When looking for the solutions to the Maxwell equations in confined 
regions (waveguides or resonators) it is convenient to use the Hertz potentials or Hertz vec- 
tors This technique considerably simplifies the calculations, specifically, removing the 
problem of choosing the gauge. It is this way that we shall follow when constructing the 
local zeta function for electromagnetic field inside a perfectly conducting wedge W a x M 1 . 

The layout of the paper is as follows. In Section II the Hertz potentials are introduced 
and the solutions to Maxwell equations inside a perfectly conducting wedge are constructed. 
These potentials are expressed in terms of two scalar functions which satisfy Dirichlet and 
Neumann boundary conditions. In Section III the local zeta function in the problem at hand 
is constructed. By a simple change of angle variable this function reduces to the analogous 
one for a massless scalar field on the manifold C^xl 1 , when j3 = 2a. Therefore the analytic 
continuation developed in Ref. || is also applicable in our case. The local zeta function 
constructed enables one to calculate the vacuum energy density of the electromagnetic field 
inside a wedge without subtraction of infinities. In Section IV transition to the global zeta 
function is carried out by introducing a cutoff nearby the cusp at the origin. Proceeding 
from this the heat kernel coefficients are calculated and the high temperature asymptotics 
of the Helmholtz free energy and of the torque of the Casimir forces are found. The wedge 
singularity gives rise to a specific high temperature behaviour ~ T 2 of the quantities under 
consideration. The obtained results are directly applicable to the free energy of a scalar 
massless field and electromagnetic field on the background of a cosmic string. In Section 
V we summarize our results. In the Appendix the local zeta functions are constructed for 
scalar massless fields inside a wedge with Dirichlet and Neumann boundary conditions. 



II. SOLUTION TO MAXWELL EQUATIONS IN TERMS OF HERTZ VECTORS 

Finding the solutions of Maxwell equations inside confined regions proves to be a rather 
difficult problem. Mainly it is due to the vector character of the electromagnetic field. A 
convenient technique for this purpose is to use the electric (IT) and magnetic (II") Hertz 
vectors P HlOU - The Lorentz gauge condition 



divA + -— ^- = 

c ot 

is satisfied identically because the Hertz vectors are connected with the potentials A and if 
in the following way 

A = -— -, y? = -divll. 

c at 

For boundary problems with cylindrical symmetry the electric (E) and magnetic (H) fields 
are expressed in terms of the Hertz vectors having only one nonzero component 
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n' = (o, o, to, n" = (o, o, v). 



(2.1) 



Here cylindrical coordinates (r, 9, z) are used with z axes directed along the intersection line 
of the wedge planes. The scalar functions U and V are the eigenfunctions of the Laplace 
operator and meet, respectively, Dirichlet and Neumann conditions on the boundary dT 



AU = u 2 U, U\ 9r = 0, (2.2) 

0. (2.3) 



AV = u 2 V, t 
on 



ar 

dV 



ar 



When considering the electromagnetic field inside a perfectly conducting dihedral of opening 
angle a (for simplicity, conducting wedge W a ) the eigenvalues of the operator (—A) are 

to 2 = oo 2 (k, A) = k 2 + A 2 , -oo < k < oo, < A < oo, (2.4) 

with k being the wave vector along the z direction. In all the functions we have dropped for 
simplicity the time-dependent factor e luJt . The spectrum ( |2.4| ) does not depend explicitly on 
the opening angle a of the wedge. This point proves to be a crucial one when attempting 
to construct a global zeta function in this problem (see Section IV). 

The Dirichlet boundary value problem ( |2.2j ) has the following normalized eigenfunctions 
(i?-modes) 

gifcz 

u Xnk (r,9,z) = —=J np {\r)s\ii{np0), 
\ hxa 



< 9 < a, p = 7r/a, n=l,2,..., < A < oo . (2.5) 



For the Neumann boundary value problem (|2.3|) the normalized eigenfunctions (if-modes) 
have the form 

v\nk{r, 9, z) = r] n0 -^-=J np (Xr) cos(np9), 
\/7ra 



— n = D 

n = 0,l,2,..., 0<A<oo, r]nv={V2' ' (2.6) 

1, n=l,2,.... 

The eigenfunctions (|2.5|) and (|2.6|) meet the orthogonality condition 

n oo poo rot -j 

dz rdr d9f* x , n , k ,(r,9,z)f Xnk (r,9,z) = 6 nn/ -6(\-\')5(k-k') (2.7) 

-oo JO JO ^ 

and the completeness condition 

°° POO POO 1 

dk \d\ f* nk (r', 9', z')f Xnk (r- 9,z) = - 6{r - r') 5{z - z') 5(9 - 9') . (2.8) 
n=0{1) J-°° Jo r 
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Here f\ n k(r,9,z) are the eigenf unctions Q2.5Q or (|2.6|) . The relations ( |2T7|) and (|2.8|) can be 
easily verified if one takes into account the formula [[l ] 



1 

A dX J q (Xr)J q {Xr') = - 6(r - r') . 

For given Hertzian vectors II' and II" the electric and magnetic fields are constructed 
by the formulas 

E = VxVxir, H = -iuo V x IT (E-modes) , 

E = iujV x II" , H = V x V x IF (//-modes) . (2.9) 



It has been proved [12| that the superposition of these modes gives the general solution to 
the Maxwell equations in the problem under consideration. 

However in the zeta function formalism there is no need to recover the electric and 
magnetic fields trough these formulas in order to construct the energy density or the energy 
momentum tensor as one proceeds when the Green's function technique is used . It is 
an essential merit of the approach considered. 



III. LOCAL ZETA FUNCTION FOR PERFECTLY CONDUCTING WEDGE 

At the beginning we remind a few general formulas that define a local zeta function. Let 
we have a differential operator L with well posed eigenvalue problem (here it is convenient 
to use the Dirac bracket notation) 

Lf„(x) = X n f n (x) or L\n >= X n \n > . (3.1) 

By making use of the completeness of the vector set { | n > } , one can represent (formally) a 
unity operator / acting in the linear space of the vectors \n > as follows 

I = ^2\n><n\. (3.2) 



71 



-1 



In view of this, we have for the inverse operator L 

L _ 1 = E \n>^ (33) 

n 

It can be easily checked with allowance for Eq. (|3.1| ). For the sth power of the inverse 
operator L^ 1 we can write 

n n 

The local zeta function ((s, x) of the operator L is a diagonal matrix element of the opera- 
tor L~ s 
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c(s,x) = 

n 

The global zeta function 



< x\n >< n\x > 
At 



n 

is obtained by integration of ((s, x) over the whole space 



as) 



x) dx . 



(3.5) 



(3.6) 



(3.7) 



Thus the local zeta function ((s,x) can be interpreted as a spatial density for the global 
zeta function ((s). 

When calculating the vacuum energy of a relativistic quantum field E by making use of 



the zeta function regularization fll^Jlifl , we have 

1 



En 



c 



(3.8) 



with the operator L being the spatial part of the relevant differential operator defining the 
field dynamics in the problem at hand. In this case A ra = u^, where u n are the classical 
characteristic frequencies of the field under study. In view of Eqs. (|3.7|) and ( |3.8|) , the 
quantity (l/2)£(— 1/2, x) can be interpreted as the vacuum energy density calculated in 
the framework of the zeta function regularization. Besides that, it is easy to show that 
the vacuum expectation value of the canonical energy-momentum tensor of non-interacting 
scalar field can be represented in the form 



< 0|T 00 (x)|0 >=\j2\ 1 n /2 f:(x)f n (x) = ~c 



1 

~2 ,X 



(3.9) 



In this case the operator L is simply —A. 

The main difficulty in obtaining the local zeta function ( |3.5| ), as well as the global one 
( [3.6D , is an analytic continuation of these formula to the region Re s < 1. 

It is remarkable that for constructing the local zeta function ((s, x) only the spectrum 
and the eigenfunctions of the relevant operator are needed. The explicit form of the regard- 
ing physical fields (for example, electromagnetic fields) is not involved in these formulas. 
Therefore we shall not use Eq. (|2.9| ) determining the electric and magnetic fields. Instead 
of that we proceed from the spectrum ( |2.4| ) in the problem at hand and the respective 
eigenfunctions (|2.5|) and ( [2.6| ). 

Substituting Eqs. (|2.4j) , (|2.5|) and (|2.6|) into the definition of the local zeta function ( p.5|) 
we obtain 



C(s,r,M = — J dk j o {k2 l\ 2)s < 3 Jg(Ar) + g J 2 np (Xr)[sm\n P e) + cos 2 (n^)] 



1 

ira 




fj (tf + A*)« 



n=l 
oo 



n=l 



(3.10) 
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Thus the local zeta function in the problem at hand is independent of the angle variable 
9 and, that is natural, of z due to the homogeneity of the system in the z-direction. For 
simplicity, we shall drop these arguments in what follows. 
By making use of the value of the integral 

dk _ -, „r(a - 1/2) 1 

- v^A 1 ' 28 v - 1 ' , Res>-, (3.11) 



,(k 2 + x 2 y v r(s) ' 2 

we cast Eq. Q3.10| ) into the form 

C(s,r) = ^^^^f:' rX 2 - 2s J 2 np (Xr)d\. (3.12) 

The prime on the summation sign means that the n = term is taken with half weight. 

Now we have to use the following integration formula (see equation 6.574.2 with fi = v 
in Ref . Ill 



r x Jl(at)dt 



a r(A)r 

-A t2, 



2i/- A + l 



Re (2i/ + 1) > Re A > 0, a > 



(3.13) 



Applying this formula we can accomplish the integration in each term in Eq. (|3.12|) for the 
respective values of the variable s restricted by the equalities 



1 < Re s < np + 



n = 0,1,2, 



(3.14) 



In these regions the integral (|3.13|) converges. Therefore the objections against using this 
formula brought up in Ref. H are obviously not applicable here. As a result, we obtain 



C(s,r) 



I r(s-l) y^,T(np + 3/2-s) 



2nar 3 2s T(s) 



n=0 



T(np + s - 1/2) ' 



(3.15) 



When deriving this equation the formula for the doubling of the gamma function argument 

m 



F(2z) = (2n)- 1 / 2 2 2z - 1 / 2 T(z)T [z + 



(3.16) 



has been used. 

By making use of the asymptotics of the ratio of two gamma functions | 16H , one can 
easily show (see below) that the series ({3 .15 ) converges in the region 



Re s > 



(3.17) 



while the first term with n = in the sum ( |3.15|) is defined in the region 



1 < Re s < - . 



(3.18) 



Thus in the problem under study, there is no representation for the local zeta function (p. 15 ) 
determining it as an analytic function of the complex variable s in a finite region. The first 
term in Eq. ( |3.15| ) is defined in the region ( |3.18| ) but the rest of the sum is defined in the 
strip 



3 3 1 

-<Re s < - + 1 . 
2 2 



(3.19) 



Cheeger has proposed to overcome this difficulty in the following way |T7| . The analytic 
continuation should be constructed separately for the first term with n = in Eq. (|3.15|) and 
for the rest sum in this equation. Further the obtained results should be added. Obviously, 
such a continuation is unique. 

Following the paper 0, we introduce the function 



G2a(s) = 



n=l 



Y{np — s + 1) 
Y{np + s) 



Re s > 1 . 



(3.20) 



For a = 7r (p = 1) we have 



T(n-s + l) 



m-s) 
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r(i-s) 

T(s) 



T(n + s) T(s) 



s(l + s) 



. . - 1 



Re s > 1, 



where F(a,f3;j;z) is the hypergeometric function. Equation 9.122.1 from Ref. [n| gives 

r(s)r(2s-2) 



F(l — s, 1; s: 1, 



r(2s - l)r(s - 1) 



Re s > 1 



Finally we obtain 



G 2 Js) 



Til 



2T s) 



Re s > 1 . 



(3.21) 



The right-hand side of this equation can naturally be considered as an analytic continuation 
of the function G2 V {s) to all over the complex plane s. On the other hand, the expression in 
the right-hand side of Eq. ( [3.21 ) for 1/2 < Re s < 1 is nothing else as the first term with 
n = in the sum in Eq. ( ^.15[ ). Thus the analytical continuation of the first term in the 
sum (|3.15| ) is trivial. 

The basic idea of the analytical continuation of the function G2 a (s) for a ^ tt proposed 
by Cheeger |T7j is also simple. It relies on the fact that the n-th term of the series (|3.20| ) 
has a polynomial asymptotics for large n 
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r ( r 7„ p It) 1) ~ WE^WW- 2 '. (3-22) 



where the functions Cj(s) are calculated exactly upon substituting in Eq. (|3.22| ) the know 
asymptotics for the gamma function From Eq. ([3.22 ) it follows immediately that the 

series in Eq. (|3.20|) determining the function G 2a (s) is convergent for 

Re s > 1 . (3.23) 

When summing over n, the terms in the asymptotics ( |3.22| ) with 

Re (2s + 2j - 1) < 1 (3.24) 

lead to the divergences. Therefore, it is natural to substitute these sums by respective values 
of the Riemann zeta function (r(s) |J 

G 2a {s) =p 1 ~ 2s P' 2j c j (s)( R {2s + 2j -1) + analytical part, (3.25) 

0<j<l-Res 

where the analytical part is the sum of the differences between a n and the asymptotic terms 
in Eq. ( 3.1 9|) picked out in the way described above. In the general case, it is impossible 
to calculate this analytical part, i.e., to express it in terms of known analytical functions. 
However, for 

11 3 h 

s = - 0,-~,-l,-^...,-~ (3.26) 



this part vanishes because for these values of s the ratio a n in Eq. ( |3.22| ) reduces to the finite 
polynomial in n. The positive integer k in Eq. ( |3.26| ) obviously depends on the number 
of asymptotic terms in Eq. ( 3.2(j| ) summed by the Riemann zeta function in the course of 
analytical continuation. Thus for values of the variable s listed in Eq. ( |3.20| ) there is an 
explicit representation of the function G2 a {s) in terms of the Riemann zeta function. For 
the boundary value s = 1 the series defining the function G 2a (s) diverges 

z — ' 1 ( np + 1 ) p L — ' n 

n=l v ' n=l 

It leads to a simple pole at the point s = 1 which acquires the analytic continuation of the 
function G 2a (s) to the region Re s < 1. 

In view of Eqs. ( |3.2(J| ) and ( |3.21 ), the zeta function ( 3.15 ) can be represented in the form 

CM = - \_J { *~ 1] [Giais - 1/2) - G 27T (s - 1/2)] = 1 3 2 / 2 " ( '" 1/2) , 

(3.28) 

where, following Ref. H, we have introduced the function 
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/ - -' (,S ^ /2) [G 2a (s) - G 2n (s) 
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[3.29) 



This function, as well as G 2a ( s ), possesses only a simple pole at the point s — 1. At the first 
sight the gamma function T(s — 1/2), entering in the definition of hais), leads to the poles for 
s = 1/2, —1/2, —3/2, .... However, it can be shown that the difference G 2a (s) — G 27T (s) 
vanishes at these points. 

The local zeta function ( |3.28| ) for electromagnetic field inside a perfectly conducting 
dihedral of opening angle a is equal to twice the corresponding zeta function for a massless 
scalar field on the space Cp x R 1 with (3 = 2a, where Cp is a cone of angle (3 ||. This 
conclusion is obvious if one takes into account the relation between corresponding boundary 
value problems (see Eq. (12) in Ref. (T^|). It is also easy to arrive at this conclusion directly 
writing out a complete set of the eigenfunction for a scalar massless field on the manifold 
Cpx 



>i. 



fxnk(r,9,z) = r) n0 - 



Akz 



=J np (Xr) Sm (np9), 

Tlfj COS 



n 



0,1,2, 



V 



2n 
J' 



(3.30) 



In this case all the eigenvalues with n = 1,2, .. . turn out to be double degenerate, each 
of them has two eigenf unctions proportional to sin(np6) or cos(np6). The local zeta func- 
tions for scalar fields inside a wedge with Dirichlet and Neumann boundary conditions are 
constructed in The Appendix. 

In order to find in the approach developed the spatial distribution of the vacuum energy 
inside the conducting wedge (see Eqs. (|3.8| ) and (|3.18| )) the value of G2a(— 1) should be 
calculated. Here we show how to do this applying the procedure of analytic continuation of 
the function G 2a (s) described above. 

From Eq. ( |3.24| ) it follows that for s = — 1 three terms with j = 0, 1, 2 in the asymptotics 
( [3.22j ) should be summed by making use of the Riemann zeta function. For s = — 1 Eq. 
( ggg ) gives 



T(np + 2) 
T(np — 1) 



[npY 



(up)' 



Comparing Eqs. (|3.31|) and fl3.22|) we infer 



co(-l) 



Cl (-1) = -1, ca(-l) 







(3.31) 



(3.32) 



According to Eq. (|3.25|) the analytic continuation of the function G 2a (s) at the point s — — 1 
has the value 



G 2a (-1) = limp 1 " 2s Vc J ( S )^( fl (2j + 2 S -i; 

s^ — 1 ' 4 



V 



j=0 

c (-l)Ci ? (-3) + c 1 (-l)p- 2 C i? (-l)+p- 4 lim c 2 (s)Cr(2s + 3) 

s^— 1 



(3.33) 



In order to find the limit in the last term in Eq. (|3.33|) , the exact form of the polynomial 
c 2 (s) should be used M 
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as well as the behaviour of the Riemann zeta function near the pole |15 

where 7 is the Euler constant. It gives 

p~ 4 lim c 2 (s)( R (2s + 3) = --^-p- A . (3.34) 
s-+-l 12U 

Substituting Eq. (|3.34j) and the values of the Riemann zeta function 

U-l) = ~ . U-3) ' 



12 ' ,jnA ' 120 
into Eq. (|3.33|) one obtains 

GW(-l) = ± (f 3 + 10 f " 7) = IS? <f 2 " + ") ■ < 3 - 35 > 
Equation (|3.21|) gives for G 27r (— 1) 

GW(-1) = 0. (3.36) 

Let us remind that upon analytic continuation of the function GWCs) the right-hand side of 
Eq. ( |3.21| ) is considered all over the plane of complex variables s. Substituting Eqs. (|3.35 ) 
and (|3.36| ) into Eq. Q3.28| ) and taking into account that 

T(-3/2) _ _2 

r(-i/2) ~ 3 ' 

we find 

<(-H = -» (p2 - i)(p2+u) (3 ' 37) 

and respectively for the density of the vacuum energy 



PW = iC(-^r) =-^^(p 2 -l)(p 2 + H), V=~ a (3.38) 



For a < 7r (p > 1) the vacuum energy density is negative, and it vanishes for a = tc (p = 1), 
when the wedge transforms into semispace. Putting in Eq. (gj3p a -> and r — > 00 m 
such a way that ra = d = constant, one arrives at the vacuum energy density between two 
parallel perfectly conducting plates 



n 2 



Pplates — 720 c/ 4 ' (3.39) 
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Multiplying this expression by the volume of space bounded by plates S d, where S is the 
area of the plate surface, we obtain the well known Casimir energy [|l9l , p0[1 



r 2 S 



E c = — ^=r • (3.40) 
G 720 c/ 3 V ; 

The vacuum energy density ( |3.38| ) posses non-integrable singularity at the point r = 0. 
Here we shall not touch the physical origin of this singularity. It is discussed in the literature 
for a long time (see, for example, Refs. 0Q,^1[]). When defining the global quantities in the 
problem at hand, we simply introduce a cutoff at the lower limit of integration over r. Thus 
for the total energy we have 



E=l rir\ HM- ^-W^ . (3.41) 



Having calculated the energy density for a conducting wedge, one can immediately re- 
cover all the components of the relevant energy-momentum tensor [Q]. In fact, taking 
into account the symmetries of the problem and the vanishing of the trace and divergence 



of Tjjj, one finds 



= ~r f(6, a) diag(l, -3, 1, 1) = (T rr , T ee , T zz , -p) (3.42) 



with some function f(8, a), vanishing when a = n. Substituting p on the right-hand side of 
Eq. (|3~42|) by Eq. fldjgp we obtain 

f{a) = T2M 2 • 

Thus all the components of the energy-momentum tensor are found. Specifically, we conclude 
that the surface density of the Casimir force F(r) acting on the wedge faces is 



(p 2 - l)(p 2 + 11) 

F(r) = -T ee (r) = [P J^J } . (3.43) 



These forces tend to diminish the opening angle a between the wedge plates. 
The torque of the Casimir forces about the origin is 

M = r Fir) r ir = ^-^^) = J_E 

J £ V ' 480vr% 2 a V ' 

Certainly it depends on the cutoff parameter e. We again recall that all the quantities found 
in our study are related to a unit length along the z direction. 



IV. HIGH TEMPERATURE ASYMPTOTICS OF THE CASIMIR EFFECT FOR 

A CONDUCTING WEDGE 

In this section we briefly consider the application of the constructed zeta function to 
investigation of the temperature behavior of the problem under consideration. In the general 
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case, to find the temperature dependence of the local characteristics proves to be a rather 
complicated problem that still waits its complete solution | 2l| . Therefore we address the 



global thermodynamical function of electromagnetic filed inside a perfectly conducting wedge 
(the Helmhotz free energy). For this purpose we need the global zeta function. However Eq. 
Q3.7D cannot be used for this purpose directly due to non-integrable singularity at r = of 
the local zeta function (|3.28 ). We again introduce a cutoff in the r integral 



CM 



i pa 

rdr / d9((s,r) 
Jo 



1 J ' S ° [G 2a (s - 1/2) - G 2n (s - 1/2)] 



2tte 1 - 2s (1 - 2s) r(s 

/ 2 q(s-l/2) 

2 y ^e 1 - 2s (l - 2s)T(s) 



(4.1) 



The appearance of the cutoff here can be formally explained in the following way. In 
the general case, the local zeta function ( |3.5| ) is, from the mathematical stand point, a 
distribution. Therefore its complete definition requires specification of the relevant set of 
test functions. In the problem at hand one can choose the test functions proportional to 
the step function 9{r — e). Certainly this consideration is a formal one, and the problem 
of constructing the global zeta function for boundary conditions in question still waits its 
solution. 

One could try to substitute the known spectrum (|2.4| ) directly in the definition of the 
global zeta function ( |3 . 7] ) . However, without introducing a cutoff the integration cannot be 
done here Really, 

c(s) - r dk r xdx - 1 v{s - i/2) r r 4 21 

c{s) -L^J (mf"v^t1 a dX - (4 - 2) 

On introducing the exponential regularization in the A-integral we obtain 

°°A 2 - 2 ^A^ [°° e- £X X 2 - 2s d\= r(3 3 r 2 ? g) - (4.3) 



o Jo 



Finally the zeta function ((s) acquires the form 



C(s) = J_ ^-i/2)r(3- 2g ) 

m 2^ T(s) • 1 ] 



This expression disagrees with Eq. ( |4.1|) . The substantial shortcoming of the formula ( |4.4| ). 
in comparison with Eq. (|4.1| ), is independence of the opening angle a. Therefore we shall 
use further the zeta function (|4.1|). 

The global zeta function (|4.1|) enables one to construct immediately the high temperature 
asymptotics of the thermodynamical functions in the problem at hand ||22|| . Here it is 
worth reminding that in our consideration we are dealing only with the spatial part of 
the differential operator governing the fluctuation dynamics of electromagnetic field (see 
Eqs. ( |2.2j ) and (|2.3|) ). Therefore the corresponding global zeta function enables one to 
construct, in a direct way, only the high temperature asymptotics of the thermodynamical 
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characteristics in the problem under study. In order to investigate complete temperature 
behavior of these characteristics a complete spectral zeta function is needed that involves 
also the summation over the Matsubara frequencies. 

The high temperature expansion for the Helmholtz free energy has the form |23| (K = 
c=l) 

F(T)~|C'(0)+ a nT 4 ~ 2n c n , T^oo. (4.5) 

n=0,l/2,l,3/2,... 

Here c n 's are numerical coefficients 

Co = ~90' Cl / 2 = -4^' Cl = "24' (4 ' 6) 

safe for c 3 / 2 and c 2 , which also depend on InT, and the heat kernel coefficients. The 

latter are determined by the zeta function through the relation 

lip ( s + n _!W)r( s ), n = 0,1/2,1,... . (4.7) 



(4tt)3/ 2 s^l-n\ 2 



As follows from Eq. ( |4.1| ) the zeta function £(s) has simple poles at the point s — 1/2 
(due to the multiplier (1 — 2s) in the denominator) and at the point s = 3/2 because the 
function I-zaiz) has such a pole at the point z — 1. Thus, only the heat kernel coefficients 
ao and a\ do not vanish in the problem under study. This is in full accord with the general 



properties of the heat kernel expansion p4| . Indeed, for flat manifolds without boundary or 
with flat boundary all the heat kernel coefficients except for ao are equal to zero. 

And what is more, one can show in the general case that the heat kernel coefficient ao is 
equal to the volume of the system under consideration. The corresponding term in the high 
temperature expansion ( |4.5|) 

°90 

is the Stefan-Boltzmann law for the free energy of electromagnetic field in the space restricted 
by the boundaries 

F it 2 2 

- = -— T 4 = --aT 4 , (4.8) 
V 90 3 V J 

where o = 7r 2 A; 4 ? /(60c 2 ^. 3 ) is the Stefan-Boltzmann constant. 

When the system has an infinite volume (just that very case is considered here), a special 
limiting procedure should be used in calculations of a (see, for example, Ref. P5| ). In the 
problem at hand this implies the introduction of a cutoff also at large values of r: e < r < R. 
In view of this the substitution of the zeta function ( f4.ip into Eq. fl4.7|) gives 

a = (R 2 - e 2 )(a - n) . (4.9) 

In the problem at hand we are dealing with two independent scalar fields. Therefore the 
value Q4.9p should be compared with doubled volume of the system under study 
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2V = 2n(R 2 - s 2 )— = (R 2 - e 2 ) a . (4.10) 

It is worth noting once more that all our formulas give the quantities per unit length along 
the axes OZ. The distinction between Eqs. ( [Op and ( |4.10| ) are obviously attributed to a 
special condition satisfied by the zeta function ([4.1|) , namely, ((s) = for a = tt. 

All the terms in the asymptotics ( |4.5D , safe for the term with n = 0, describe the deviation 
from the Stefan-Boltzmann law. In the problem in question the coefficient a\ does not vanish 
due to the discontinuity of the boundary at the point r = (wedge or cone singularity). 
There is a voluminous literature concerning with such a singularity, we mention here only 
a few papers [|6l, |T7t^6H 29|| where further references can be found. 

Taking into account the behaviour of the gamma function at the origin r(s) ~ s — > 

we deduce from Eq. Q4.1D 

= M-V 2 ) . (4.ii) 

The explicit formulas derived in Ref. give 



,,,-,2) = -^ = ^. 1112; 



Hence 



CO) = (4.13) 



Substituting the zeta function ( f4.1p in Eq. ( [4.7] ) we find 



ai lim (s - l/2)C(s)r(s) = -^4^- ( 4 - 14 ) 



(4tt)3/2 s _ > i/ 2 V 4^ 

In order to calculate (0) we take, in accord with relation ( |3.24| ), two terms in Eq. ( |3.25| ) 
with j = 0,1. It gives 

G 2a (0) = \imp 1 - 2s [co(s)CR(2s-l)-p- 2 c 1 (s)C R (2s + l)}. (4.15) 

s— >0 



The coefficients cq(s) and ci(s) are explicitly written, for example, in Ref. |TJ 



cb(s) = 1, Cl (s) = |(s - l/2)(s - 1) . (4.16) 

With allowance for the behaviour of the Riemann zeta function close by pole 

( R (l + 2s) ~ -L + 7 + ..., s ^0 
2s 



we obtain 



GWOH^p- 1 -?). (4.17) 
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Form Eq. ( |CTD it follows that 



Finally we have 



and 



GW(0)=0. (4.18) 



1 , _ lN 1 fit a 



/ 2a (0) = -(p-p- 1 ) = - 

D 6 \Q 71 



ai = -27rJ 2a ) = . (4.19) 

3 Va 7T/ 

Now we are in position to write the high temperature corrections to the Stefan-Boltzmann 
law inside a conducting dihedral angle 

F cor {T) ~ -^«@) + ^-^T 2 , T^oc. (4.20) 
corV ; 32tt% 72aftc ' V ; 

Here the constants h and c are explicitly restored in order to separate a pure classical 
contribution (the first term on the right hand side of Eq. (|4.20|) ) and the quantum correction 
(the second term in Eq. ( 4.20|) ). 



The authors of Ref. |7j argue that the high temperature asymptotics of the free energy 
of electromagnetic field inside a perfectly conducting wedge should be the same as that for 
smooth conducting surfaces with nonzero curvature 

F(T) ~ CT\n(T/Q), T -> oo, (4.21) 

with C and Q being some const ant sQ. The specific term proportional to T 2 in the asymptotics 
( [4.21]) was not brought out. 

Making use of Eq. ( |4.20 ) one can try to estimate, at least in kind, the behaviour of the 



torque of the Casimir forces at high temperature. For this purpose we substitute the total 
energy E in Eq. (|3.44|) by the asymptotic expression for the free energy in the problem under 
study 

_ 3( R (3)T 7i 2 -a 2 T 2 „ , . 

Af~ V 5-1^, T^oo. 4.22 

32vr 3 ae 24a 2 he' K J 

The book [0 is a sole source where the high temperature behaviour of the torque M has 
been considered. The authors have used the method of images for constructing the Green's 
functions of electromagnetic field inside the dihedral, the photon energy being cut in ultra- 
violet region. In terms of the notations used in the present paper, the high temperature 
asymptotics of M found in Ref. Jl| reads 

M ~ — , T -> oo. (4.23) 

12a V V 1 



x The examples of such asymptotics can be found in Ref. [^] 
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Here we have also to ascertain disagreement with our calculations. 

In the preceding section it was shown that the zeta function of electromagnetic field 
inside a perfectly conducting dihedral of opening angle a is equal to double zeta function 
for a massless scalar field on the manifold Cp x M 1 with 2a = /3, where Cp is a cone of angle 
(3. Therefore all the results obtained in our consideration are directly applicable to the 
massless scalar field on the manifold Cp x K 1 , in particular to the cosmic string background. 
Really, the general solution to the Maxwell equations in the presence of a cosmic string 
are expressed in terms of two scalar massless fields that satisfy the periodicity conditions 
in angle variable with a period j3 = In — <f), where <fi is the deficiency angle of the metric 
due to the cosmic string, <p = 87tGfi. Most simply one can come to this conclusion in 
the following way. Let us assume that the string has a finite transverse size, r str , and its 
surface is perfectly conducting. Obviously, the reasonings of the Section II apply to this 
boundary configuration, two scalar functions U and V being periodic in angle variable 9 
and meeting the Dirichlet and Neumann boundary conditions at r = r str . Letting r str 
tend to zero, we substitute the Dirichlet and Neumann boundary conditions at r = r str 
by requirement that the relevant solutions to be bounded everywhere including the origin: 
< r < oo. The periodicity condition in angle variable 9 apparently remains. From here we 
infer that the local zeta function for electromagnetic field in the presence of cosmic string is 
obtained from the relevant zeta function for the conducting wedge (|3.28| ) by the substitution 
2a = 2n — 87rGfi. Specifically, for the vacuum energy density in this case we have again Eq. 
( p.39| ) with the same substitution. The same result was obtained in Ref. |3(J by making use 



of the Greens function method. Further, the high temperature corrections to the Stefan- 
Boltzmann law for the free energy of electromagnetic field on the background of a cosmic 
string are given by the formula 

Fm *™&+*L + *=4, T^oo. (4.24) 
v ; 32vn% 72/ic2 2vt -6' V ; 



V. CONCLUSION 

We have shown that the technique of the local zeta function enables one to carry out, in 
a consistent way, the calculation of the vacuum energy density of electromagnetic field inside 
a conducting wedge without dealing with obvious divergences. This regularization method 
leads to a finite expression for the energy density under study without any subtractions. 
Our consideration is a natural addition to studies of the Casimir effect for a wedge accom- 
plished in recent papers [^|J by making use of Green's function method. Employment of 
the Hertz potentials for constructing the general solution of the Maxwell equations results in 
a considerable simplification of the calculations of the local zeta function in the problem at 
hand. Transition to the global zeta function has been carried out by introducing a cutoff for 
small r. On this basis the nonzero heat kernel coefficients do, Or, and the zeta determinant 
were calculated. Remarkably that the coefficient at does not depend on the cutoff and it is 
completely determined by the wedge singularity of the boundary. Proceeding from this the 
high temperature asymptotics of the Helmholtz free energy and of the torque of the Casimir 
forces are found. The wedge singularity gives rise to a specific high temperature behaviour 
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~ T 2 of the quantities under consideration. The found results are directly applicable to the 
free energy of electromagnetic field on the background of a cosmic string. 
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APPENDIX: LOCAL ZETA FUNCTIONS FOR SCALAR MASSLESS FIELDS 
INSIDE DIHEDRAL ANGLE WITH DIRICHLET AND NEUMANN BOUNDARY 

CONDITIONS 



First we consider the Dirichlet boundary conditions with the eigenfunctions 
local zeta function is given by 



The 



(d(s,t, 



OO POO 



— I (W 



oo JO 



XdX 



(k 2 + s 



2\s 



XX(Ar) sin 2 (n P 9) 



r(*-l/2) 
2^aT(s) J 



n=l 
oo 



dXX 2 - 2s Y,Jn P W[l-cos(2np)}. 



(Al) 



n=l 



By mens of the formula 8.531.3 from [15 ] 

oo 

J (zsma) = Jl (|) + 2^ J 2 k (|) cos{2ka) 



k=l 



the dependence on the angle 9 in Eq. (Kl) can be stored in a single term 



( D (s,r,9) 



r(a - 1/2) 
2y/ifaT(s) Jo 



dXX 



2-2s 



oo 

- J 2 (Xr) + Jj,(Ar) - -J (2Ar sinp^) 



n=l 



2 US ' } A^aT Jo 



dXX 2 - 2s Jo(2Xrsm( P 9)), 



(A2) 



where ((s, r) is the local zeta function for electromagnetic field in a conducting wedge defined 
in Eq. ( |3.10| ). By making use of the integration formula 6.561.14 from [|15| 

n T (~ + ~z/ + 1 
x fl JJax)dx = 2^a-^ x — )| f f-f -Rei/-l<Rew<-, a>0 

r (l + ^-^) 2 



we cast Eq. (|A2|) into the final form 
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Us, r, 9) = Ids, r) + - 1 F(3 ^ s) . (A3) 

The last term in Eq. ( |A3| ) is first defined in the strip 

< Re s < |, (A4) 

determined by the conditions of applicability of Eq. (|A3|). As usual under the analytic 
continuation, we assume that outside this strip the same equation holds. 

Proceeding along similar lines we obtain the local zeta function for the Neumann bound- 
ary conditions 

, \ , , * 1 r(3/2-s) /Ar , 

C D (s,r,9) = -C(s,r) - — — . 2g 1 ' ' . A5 

In particular the vacuum energy densities for massless scalar fields subjected to the Dirichlet 
and Neumann boundary conditions on the wage sites are 

p D (r,9)=p_(r,6), p N (r,9) = p + (r,6), (A6) 

where 

p±(r,9) = \p(r)± / j- - , (A7) 

2 167rar 4 sin yp9) 

and p(r) is the vacuum energy density of electromagnetic field, defined in Eq. ( |3.38| ). 

Thus the vacuum energy densities derived have nonintegrable singularities near by the 
wedge sites. It is a typical behaviour of spatial distribution of the Casimir energy (see, for 
example, Ref. |j3~I]). 

The local zeta function technique, consider in the present paper, cannot directly be 
applied to the scalar field with conformal coupling. The point is that in this case the 
relation (B.9I) does not hold. 
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